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This paper reports on further results concerning nonlinear equations 
of the form F(x) + Ax = B, in which F(-) is a "diagonal nonlinear 
mapping" of real Euclidean n-space E" into itself, A is a real n X n 
matrix, and B is an element of E n . Such equations play a central role in 
the dc analysis of transistor networks, the computation of the transient 
response of transistor networks, and the numerical solution of certain 
nonlinear partial-differential equations. 

Here a nonuniqueness result, which focuses attention on a simple special 
property of transistor-type non-linearities, is proved; this result shows that 
under certain conditions the equation F(x) + Ax = B has at least two 
solutions for some B (E E". The result proves that some earlier conditions 
for the existence of a unique solution cannot be improved by taking into 
account more information concerning the nonlinear ities, and therefore 
makes more clear that the set of matrices denoted in earlier work by P 
plays a very basic role in the theory of nonlinear transistor networks. In 
addition, some material concerned with the convergence of algorithms for 
computing the solution of the equation F(x) + Ax = B is presented, and 
some theorems are proved which provide more of a theoretical basis for the 
efficient computation of the transient response of transistor networks. In 
partindar, the following proposition is proved. 

If the dc equations of a certain general type of transistor network possess 
at most one solution for all B £ E" for "the original set of a's as well as 
for an arbitrary set of not-larger as", then the nonlinear equations en- 
countered at each time step in the use of certain implicit numerical inte- 
gration algorithms possess a unique solution for all values of the step size, 
and hence then for all step-size values it is possible to carry out the algo- 
rithms. 



* The material of this paper was presented at the Advanced Study Institute on 
Network Theory (sponsored bv the N.A.T.O.; Knokke, Belgium; September 1-12, 
1969). 
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I. INTRODUCTION AND DISCUSSION OF RESULTS 

References 1 and 2 present some results concerning the equation 

F(x) + Ax - B, (1) 

in which, with n an arbitrary positive integer, A is a real n X n matrix, 
B is an element of real Euclidean n-space E n , and F(-) is a mapping of 
E n into E n denned by the condition that* for all x = fa , x 2 , • • • , x n ) ir E 

E n , 

F(x) = tfjW,/^), ••• ,/ n (.Or (2) 

with each /,(•) a strictly monotone increasing mapping of E 1 into itself. 
Equation (1) plays a central role in the dc analysis of transistor net- 
works,** the transient analysis of transistor networks (see Section 1.4), 
and the numerical solution of certain nonlinear partial differential 
equations. 

In Ref . 1 it is proved that there exists a unique solution x of equation 
(1) for each strictly monotone increasing mapping F(-) of E" onto E n 
(that is, for each set of strictly monotone increasing mappings /,(•) of 
E 1 onto itself) and each B E E n if and only if A is a member of the set 
P of real nX n matrices with all principal minors nonnegative. It is also 
proved in Ref. 1 that equation (1) possesses a unique solution x for 
each continuous monotone nondecreasing mapping F(-) of E" into E n 
(that is, for each set of continuous monotone nondecreasing mappings 
of E l into E 1 ) and each B E E n if A belongs to the set P of all real 
n X n matrices with all principal minors positive 1 '. A direct modification 
of the existence proof given in Ref. 1, as indicated in Ref. 2, shows that 
equation (1) possesses a unique solution for each strictly monotone 
increasing mapping F(- ) of E n onto (a t , ft) X (a 2 , /3 2 ) X ■ • • X (a„ , /3„) 
with each a, and /3, elements of the extended real line (real line) such 
that a,- < 0i and each B E E n if (and only if) A G Po and det A ^ 0. 
Some network theoretic implications of these and related results are 
discussed in Refs. 1 and 2, where the matter of determining whether or 
not A E Po or A E P is considered in some detail. 



* Throughout the paper the superscript tr denotes transpose. 

** See Ref. 1 for a derivation of the equation within the context of the transistor 
dc-analysis problem. 

t There are some interesting applications of this result in the study of numerical 
methods for solving certain nonlinear partial-differential equations, in which A 
has nonpositive off-diagonal terms and is irreducibly diagonally dominant. 3 

* The numbers a ( and /3 { are members of the extended real line if — oo ^ on 
S oo and — oo < fl ( ^ oo. 
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This paper presents a proof of a nonimiqueness result. The proof 
focuses attention on a simple special property of transistor-type non- 
linearities. The result shows that under certain conditions equation (1) 
has at least two solutions for some B £ E n . In addition, the paper 
presents some material concerned with the convergence of algorithms 
for computing the solution of equation (1) and proves some theorems 
which provide more of a theoretical basis for the efficient computation 
of the transient response of transistor networks. The remaining portion 
of Section I is concerned with a detailed discussion of the results and 
their significance. 

1.1 An Application of the Nonuniqueness Theorem 

The standard Ebers-Moll transistor model, which is widely used, 
gives rise to functions /,(•) which, while continuous and strictly mono- 
tone increasing, are mappings of E x onto open semi-infinite intervals. 
For such /,■(•)> the results stated above assert that the equation (1) 
possesses at most one solution x for each B £ E n if A £ P ; and if 
A £ P and det A =^ 0, then equation (1) possesses a solution for each 
B £ E n . Since, as indicated in Ref. 1, A = T~*G with T a nonsingular 
matrix which takes into account the forward and reverse transistor a's, 
and G is the short circuit conductance matrix of the linear portion of 
the network, the condition that det A not vanish is equivalent to the 
rather weak assumption that the linear portion of the network possess 
an open-circuit resistance matrix. 

It is natural to ask whether the use of more-detailed information 
concerning the nonlinearities of the transistor model would enable us 
to make assertions concerning the existence of a unique solution of 
equation (1) for all B £ E n under weaker assumptions on A. In particu- 
lar, can the condition that A belong to P be relaxed? The first result 
proved in this paper, Theorem 1 of Section II, shows that if the /,(•) are 
exponential nonlinearities of the type associated with the Ebers-Moll 
model, then the condition that A belong to P cannot be replaced by 
a weaker condition. More explicitly, in Section II a set T Q of mappings 
of E" into E" is defined, and T contains all of the mappings F(-) that 
correspond to Ebers-Moll type /,•(• )'s. It is proved there that if A (£ Bo , 
then for any F(-) £ So , there is a B £ E n such that equation (1) 
possesses at least two solutions. In fact, it is proved that if A (£ P 
and if 5 is an arbitrary positive number, then for any F(-) E 3o i there 
is a B £ E" such that equation (1) possesses two solutions such that 
the distance in E n between the two solutions is 5. 

Thus Theorem 1 together with the earlier results mentioned above 
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concerning existence of solutions show that the set of matrices P plays 
a quite fundamental role in the theory of nonlinear transistor networks. 

1.2 An Algorithm for Computing the Solution of Equation (1) 

Several results which assert that A £ P under certain conditions on 
the transistor as and the short-circuit conductance matrix of the linear 
portion of the network are proved in Refs. 1 and 2. In particular, 
Ref. 1 proves that A £ P, and hence that A G Po , if A = P~*Q with P 
and Q real n X n matrices such that for all ; = 1,2, • ■ • , n 

Pa > Z) \Va I and Qa > £ I 9u \* 

Theorem 2 of Section II shows that a relatively simple and entirely 
constructive algorithm can be used to generate a sequence x m , x (1 , 
of elements of E n that converges to the unique solution of (1) if A = 
P~ l Q with P and Q as denned above and each /,(•) is a continuous (but 
not necessarily differentiable) monotone nondecreasing mapping of 
E 1 into E 1 ** 

1.3 Palais' Theorem, Existence of Solutions of Equation (1), and Algo- 
rithms for Computing the Solution of Equation (1) 

Reference 1 gives two existence proofs concerning equation (1). One 
proof, the more basic of the two, is based on first principles and em- 
ploys an inductive argument in which, with k an arbitrary positive 
integer less than n, the existence proposition is assumed to be true with 
n replaced by k and it is proved that then the proposition is true with n 
replaced by (A; -f 1). The second proof uses a theorem of R. S. Palais 
and requires that the /,-(•) be continuously differentiable throughout E 1 . 
More explicitly, Palais' theorem 1 " asserts that if R(-) is a continuously 
differentiable mapping of E n into itself with values R(q) for q E E n , 
then jR(-) is a diffeomorphism : of E n onto itself if and only if 

(i) det J Q 7* for all q £ E tt , in which J Q is the Jacobian matrix of 
R(-) with respect to q, and 

(it) \\R(q) || -> co as || q || -* •« 

* It is proved also that A 6 P if A = P~ l Q with p„- > X.vy | p,-,- land?,,- ^ 

** A related result given in Ref. 4 is not directly applicable here because of assump- 
tions made in Ref. 4 concerning the existence and boundedness of a certain Jacobian 
matrix. 

t See Ref. 5 and the appendix of Ref. 6. 

* A diffeomorphism of E n onto itself is a continuously differentiable mapping of 
E n into E" which possesses a continuously differentiable inverse. 

tt Here || • || denotes any norm on E". 
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And the second proof of Ref . 1 shows that, with 

R(q) = F(q) + Aq 

for all q E E n , the two conditions (i) and (ii) are met when A E -Po and 
each /,(•) is a continuously differentiable strictly-monotone-increasing 
function which maps E 1 onto E l and whose slope is positive through- 
out E 1 .* 

There are some problems which arise in connection with, for example, 
the numerical solution of certain nonlinear partial-differential equations** 
in which one encounters an equation of the form (1) with A ^ P and 
det A ?* 0, but with functions /,(•) which, while continuously differ- 
entiable, are monotone nondecreasing (rather than strictly monotone 
increasing) mappings of E 1 into E l . We can prove that even in such 
cases equation (1) possesses a unique solution for each B E E n as follows. 
Here the Jacobian matrix of F{q)-\-Aq exists and is of the form D(q)-\-A 
in which D (q) is a diagonal matrix with nonnegative diagonal elements. 
Since A E P and det A j* 0, we have 2 det [D(q) + A] j* for all 
q E E n . An immediate application of Theorem 3 of Section II shows 
that || F(q) -f Aq \\ — > «. as || q || — > «.* Therefore, by Palais' theorem, 
F(x) + Ax = B possesses a unique solution for each B. 

Theorem 3 is of use not only in connection with the proof given in 
the preceding paragraph; it also plays a key role in showing that there 
is an algorithm which generates a sequence of elements of E n x ,x , • • • 
that converges to the unique solution of F(x) + Ax = B whenever each 
/<( ■ ) is twice continuously differentiable on E 1 and the conditions on A 
and F(-) of the preceding paragraph are satisfied. 1 

More generally, if R(-) is any twice-differentiable mapping of E n 
into itself such that conditions (i) and (ii) of Palais' theorem are 
satisfied, then, with R~ l (-) the continuously-differentiable inverse of 
R(-), x = R~\B) satisfies R(x) = 9 in which is the zero element of E n , 
and there are steepest decent as well as Newton-type algorithms each of 



* The reasons that two proofs were presented in Ref. 1, with the second proof 
a proof of a somewhat weaker result, are that the arguments needed for the appli- 
cation of Palais' theorem had already been developed in Ref. 1 and used for 
other purposes there, and it was felt desirable to indicate an alternative approach 
to essentially the same problem. 

** The writer is indebted to J. McKenna and E. Wasserstrom for bringing this 
fact to his attention. 

f More explicitly. Theorem 3 shows that there is a vector C G E" such that 
|| p(q) _|_ ^io _(_ c'\\ — * oo as || g II — » °°, which is equivalent to the statement 
concerning |j F(q) + Aq || made above. 

* The differentiability assumption here is introduced as a matter of convenience, 
and is certainly satisfied when the /,(•) are Ebers-Moll exponential-type non- 
linearities. 
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which generates a sequence in E" that converges to x. To show this, let 7 
j(y) = R(y) || 2 for all y £ E n in which || • || denotes the usual Euclidean 
norm (that is, the square-root of the sum of squares) . Since condition (j) 
of Palais' theorem is satisfied, the gradient V/ of /(•) satisfies (V/)(i/) 
5^ unless f(y) = 0,* and since condition (ii) of Palais' theorem is satis- 
fied the set S = [y £ E n : j(y) ^ f(x m } is bounded for any x M £ E". 
Therefore we may appeal to, for example, the theorem of page 43 of Ref. 
7 according to which for any .t (0> £ E n , for any member of a certain class 
of mappings <£>(•) of S into E", and for suitably chosen constants y n , 
7i , • • • , the sequence .r (0) , .r (1) , • • • defined by 

x (* + « = x (« + TtVJ ( x <«) for all A- ^ 

belongs to S and is such that || R(x m \\ — * as k — > w. However, since 
/^~'(-) exists and is continuous, 1 it follows from 

x (k) = R- 1 [R(x w )] for all k^O 

and the fact that R(x (k) ) — > 6 as A; — » w, that limjt..,,, .r a) exists and 

lima; <w = iT 1 ^), 

A-«oo 

which means that x = limj i _ 00 a; ( * } . 

1.4 Transient Response of Transistor-Diode Networks and Implicit Nu- 
merical-Integration Formulas 

At this point we briefly consider some aspects of the manner in which 
the previous material bears on the important problem of providing 
more of a theoretical basis for numerically integrating the ordinary 
differential equations which govern the transient response of nonlinear 
transistor networks. Although we consider explicitly only networks con- 
taining transistors, diodes, and resistors, the material to be presented 
can be extended to take into account other types of elements as well. 
In addition, we shall focus attention on the use of linear multipoint 
integration formulas of closed (that is, of implicit) type, since such 



* Here we have need the fact that ( V/)(2/) = 2J v ' r R(y) for all y G #". 7 
t By Palais' theorem R ( • ) is a diffeornorphism of E" onto itself. 

* The material of the second part of Section 1.3 was motivated by previous 
recent work of the writer's colleague A. Gersho who made the observation that the 
convergence of an algorithm for the solution of equation (1) could be shown by 
combining results of Ref. 1 with the approaches described by Goldstein. 7 (See the 
November 1969 B.S.TJ. Brief by A. Gersho.) 
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formulas are of considerable use in connection with the typically "stiff 
systems" of differential equations encountered. 

A very large class of networks containing resistors, transistors, and 
diodes modeled in a standard manner is governed by the equation 8 

^ + TF[C~ l (u)] + (/ + GRY'G€r\u) = B(t), t^O (3) 
at 

where, assuming that there are q diodes and /) transistors, 

(i) T = I q © I\ © T 2 © • • ■ © T p , the direct sum of the identity 
matrix of order q and p 2 X 2 matrices T k in which 

r. = l — 

with < a ( r k) < 1 and < a) k) < 1 for k = 1, 2, ■ ■ ■ , p. 

(ii) R = R n © #, © R 2 © • • • © R P , the direct sum of a diagonal 
matrix R = duxg (r, , r 2 , • • • , r a ) with r k ^ for k = 1,2, ■ • • , q 
and p2 X 2 matrices 72 A in which for all k = 1, 2, • • • , p 

A-K +f{M ^ 1 

with r"' ^ 0, rj" ^ 0, and rj w ^ 0. (The matrix # takes into account 
the presence of bulk resistance in scries with the diodes and the emitter, 
base, and collector leads of the transistors.) 

(Hi) G is the short-circuit conductance matrix associated with the 
resistors of the network. (It does not take into account the bulk re- 
sistances of the semiconductor devices.) 

(iv) F(-) is a mapping of E (2p+9) into E (2p+Q) defined by the condition 
that 

F(x) = [/,(.r,),/ 2 (.r 2 ), ••• , / 2 „ +a (.r 2/)+ ,)]' r 

for all x £ E (2p+Q) with each /,(•) a continuously differentiable strictly- 
monotone increasing mapping of E l into E 1 . 

(v) C _1 (-) is the inverse of the mapping C(-), of E (2p+q) into itself, 
defined by 

C(x) = diag (d , c 2 , • • • , c 2p+Q )x + diag (t x ,t 2 , • ■ ■ , t 2p+q )F(x) 

for all x G E {2p+q) with each c, and each r,- a positive constant. 

(vi) B(t) is a (2p + </)-vector which takes into account the voltage 
and current generators present in the network, and 
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(vti) u is related to v the vector of junction voltages of the semi- 
conductor devices through C(v) = u for all v £ j g <2p+<i > 
Equation (3) is equivalent to* 

u + f(u, t) = (2p+a) , t^O (4) 

in which of course 

f(u, t) = TF[C-\u)] + (I + GRY'GC-'tu) - B(t) (5) 

and 0( 2p+ a) is the zero vector of order (2p + <?)• 

It is well known that certain specializations of the general multi- 
point formula 9,10 

r r 

2/n+l = 2 VkVn-h +'l E M»-t (6) 

k=0 k=—l 

in which 

ft.-* = -flVn-k , (n - k)h] (7) 

can be used as a basis for computing the solution of equation (4). 
Here h, a positive number, is the step size, the a k and the b k are real 
numbers, and of course y n is the approximation to u(nh) for n ^ 1. 
In the literature dealing with formulas of the type (6) in connection 
with systems of equations of the type (4), information concerning the 
location of the eigenvalues of the Jacobian matrix J u of f(u, t) with 
respect to u plays an important role in determining whether or not a 
given formula will be (in some suitable sense) stable. In particular, an 
assumption often made is that all of the eigenvalues of J u lie in the 
strict right-half plane for all t ^ and all u. For f(u, t) given by equa- 
tion (5), we have 

/.-rdh« { ffl*M J 

\c t + T,/,[0,(W,)]J 

+ (7 + aB )- g duH t { e< + ^ [g)(U()] } (8) 

in which for j = 1,2, • • • , (2p + q) ^.(w,) is the j tb component of C -1 (w). 
Thus here «/„ is a matrix of the form 

TD, + (J + GR)- l GD 2 (9) 

where Dx and Z) 2 are diagonal matrices with positive diagonal elements. 



* Ref. 8 shows that if B( • ) is a continuous mapping of [0, «> ) into E^p + "\ then for 
any initial condition w (0) £ £0p+«) there exists a unique continuous (2p + q)- vector- 
valued function u(-) such that «(0) = w (0) and (3) is satisfied for all t > 0. 
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A simple result concerning equation (9), Theorem 4 of Section II, 
asserts that if there exists a diagonal matrix D with positive diagonal 
elements such that* 

(i) DT is strongly column-sum dominant, and 
(ii) D(I + GR)~*G is weakly column-sum dominant, 

then for all diagonal matrices D x and D 2 with positive diagonal ele- 
ments, all eigenvalues of (9) lie in the strict right-half plane. This con- 
dition on T, G, and R is often satisfied. 1 " 

The subclass of numerical integration formulas (6) defined by the 
condition that b- r > are of considerable use 11 ' 12,13 in applications 
involving the typically "stiff systems" of differential equations en- 
countered in the analysis of nonlinear transistor networks. With &_i > 0, 
y n+l is defined implicitly through 

r r 

y n+l + hb-if(y n +i , (n + l)/i) = X) a *#»-* + h X) &*&-* 

*=0 *-0 

in which the right side depends on y n - k only for A; £ {0, 1, 2, • • • , r\, 
and for f(u, t) given by equation (5), we have 

fr+i + »-i{r^[(T , (y. +1 )] + (J + GRy'GC-^y^)) = <7„ (10) 
in which 

q n = E a k y n . k + h £ 6 A .y n _ fc + *6_,B[(n + 1)*]. 

Obviously, the numerical integration formula (10) makes sense only 
if there exists for each n a y B+1 £ .E <2p+a> such that equation (10) is 
satisfied. 

Let .T n+1 = C~ l (y n+1 ) for each w. Then equation (10) possesses a 
unique solution y n+1 if and only if there exists a unique x n+1 £ i? 
such that 

C(o; n+1 ) + to_ 1 [r/?(a; B+l ) + (/ + GAr^i] = g„ . (11) 
Since C(x n+i ) = r.r n+ , + rF(x n+i ), in which 

c = diag (c, , c 2 , • • • , c 2p+a ) 
and 

t = diag (t, , r 2 , ••• , t 2p+ „), 



* The terms "strongly -colli rrm-sum dominant" and "weakly -column-sum domi- 
nant " are reasonably standard. However they are defined in Section II. 
t See Ref. 8 for examples. 
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equation (11) is equivalent to 

[r + hb^T]F(x n+1 ) + [c + *MJ + ffi^-'flteM., = <z„ ■ (12) 

The matrices r and c are both diagonal with positive diagonal ele- 
ments. Thus it is clear that for all positive /* 

det [r + hb-xT] ?± 

and 

det [c + Kb-td + GR)~ l G\ * 0.* 

For all sufficiently small positive h 

[r + hb^T]- l [c + hb^il + GRY'G] £ Po - + 

Consequently ft for all sufficiently small /; > 0, equation (12) possesses 
a unique solution for each q„ . x However, our interest in equation (12) 
is primarily in connection with "large-/*" algorithms. 

Suppose that det G ^ and that T~*G £ -Po for all possible combi- 
nations of a r and a, (0 < a T < 1,0 < a s < 1) for each transistor (see 
Ref. 1 for examples). Then, according to Theorem 6 of Section II, 
for any particular T and R 

[r + hb^TT l [e + AM/ + GRy'G] £ Po 

for all h > 0, and hence equation (10) possesses a unique solution 
y n+1 for all positive values of /*. 

An important and general proposition concerning (10) is as follows. 
Suppose that 

T-'[(J + GRY l G] £ Po (13) 

and that condition (13) is satisfied whenever a[ k) and a\ k) are replaced 
with positive constants 8 l r k) and 8 ( f k) , respectively, such that 5'*' ^ a {k) 
and 8 lk) ^ a {k) for k = 1, 2, • • • , p. In other words, assuming that P(-) 
is as denned in this section and that F(-) £ Jo 2 "*"' ( see Definition 1 of 
Section 2.1), suppose that the dc equation 

F(x) + T~ l [(I + Gfl) _1 (?].T = B 

possesses at most one solution x for each B £ E 12 "*^ for "the original 
set of a's as well as for an arbitrary set of not-larger as." Then an 

* Here we have used the fact that (J + GR)~ l G is positive semidefinite. 
+ See Section 1.2. 

tt See Section 1.3. 

* Alternatively, this conclusion could have been obtained by applying the con- 
traction-mapping fixed-point principle to (10), in view of the fact that each of the 
elements of J u is bounded on u £ £<*+«) and t £ [0, oo ). 
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immediate application of Theorem 5 of Section II shows that 
[t + hb-xTT^c + hb-M + GRy'G] G Po 

for all h > 0, and hence that equation (10) possesses a unique solution 
i/ n+1 for all A > and all q n G E l2v+Q) . 

II. THEOREMS, PROOFS, AND SOME DISCUSSION 

Throughout this section, 

(i) n is an arbitrary positive integer, 

(ii) P Q denotes the set of all real n X n matrices M such that all 
principal minors of M are nonnegative, 

(Hi) real Euclidean ??-space is denoted by FT, and 6 is the zero 
element of E", 

(iv) v tr denotes the transpose of the row vector v = (v l ,v 2 , • • ■ ,v„), 
(v) || v || denotes (2/-i W for aU v G E n , 

(vi) if D is a real diagonal matrix, then D > (D ^ 0) means that 
the diagonal elements of D are positive (nonnegative) , 

(pit) I„ denotes the identity matrix of order q, and / denotes the 
identity matrix of order determined by the context in which the symbol 
is used, and 

(to) we shall say that a real n X n matrix M is strongly (weakly) 
column-sum dominant if and only if for j = 1,2, • • • , n 

ma > (^) X I m a I- 

2.1 Definition 1 

For each positive integer n, let $o denote that collection of mappings 
of E" into itself denned by: F £ 3^ if and only if there exist for j = 
1, 2, • • • , n, continuous functions /,(•) mapping E l into E 1 such that 
for each x = (.r, ,x 2 , • ■ • , .r„)' r G E n , F(x) = [fifa), / 2 (x 2 ), • • • , /„(x n )]' r , 
and 

(i) inf [/,-(« + /3) - /,-(«)] = 

o«(-c».oo) 

(ii) sup [/,-(« + /3) - /,-(«)] = + co 

o< ( — oo.co) 

for all |8 > and all j = 1, 2, • • • , n. 

2.2 Theorem 1 

Let F £ $1 , let A be a real n X n matrix such that A (£ P , and let 
b be a positive constant. Then there exist B G FT, x G FT, and y G FT 
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such that 

(i) F(x) + Ax = B, 

iii) F{y) + Ay = B, 

and 

{Hi) II x - y || = 8. 

2.3 Proof of Theorem 1 

Since A ££ P , there exists 2 a real diagonal matrix D > such that 
det (D + A) = 0. Thus there exists a x* £ _#" such that || x* || = 5 
and (D + ^L)x* = B. 

Since -P £ 5q , there exists &x (E. E" such that 

/*(*#) -/#(* -**> = «1* 

for all ;' = 1,2, • • • , n in which d { is the j* diagonal element of D. Let 

B = F(x) + At, 
and let y = x - x*. Then A(x - y) = Ax* = -Dx*, and 
F(x) - F{y) + A(x - y) = 6. U 

2.4 Remarks Concerning Theorem 1 

If, as in the case of standard transistor models, 
/,.(*,) = e*« - 1 

or 

/,(*,) = 1 - e-*« 

with X, > 0, we have, respectively, 

or 

/,-(« + /3) - /,(«) = 9~ ltm (1 - e- x ' fl ) 

and it is clear that for either type of function conditions (i) and (ii) 
of Definition 1 are satisfied. 

In Ref. 1 it is proved that if F(- ) £ 3TZ the set of all F( ■ ) of the form 
(2) with each /,(•) a strictly monotone increasing mapping of E 1 into 
E 1 , and if A £ P , then equation (1) possesses at most one solution. 
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Thus, using Theorem 1, we see that for each F(-) £ 3TC (~\ ffJJ there 
exists at most one solution of F(x) + Ax = B for each B £ W if and 
only if A £ P . Similarly, with 3TCo the set of all F(-) of the form (2) 
with each /,(•) a strictly monotone increasing mapping of E 1 onto 
(a,- , fit) with each a, and (J { such that — °o ^ a, < /Si ^ <» , if P(-) £: 
9TCo O Jo and det A j£ 0, then there exists a unique solution x of 
P(x) + Ax = B for each 5 £ # n if and only if A G P • (The "if" part 
of this statement is proved in Ref. 2.) A parallel development can be 
carried out for equations of the form AF(x) + x = B with A a real 
n X n matrix, F( • ) £ 9IZo H 9^, and B £ #". More explicitly, we can 
prove that if F( • ) £ 9TC O ffo > then there exists a unique solution x of 
AF(x) + x = B for each P. £ E n if and only if A £ P . 

There may be a temptation to conjecture that whenever F(-) £ 9TC C\ 
5o and A £j: P then the equation F(x) -{-Ax = B does not possess a 
solution for some B £ E". The conjecture is false. In fact, with n = 2, 
A(xi) = e", / 2 (x 2 ) = e x ", and 



A = 



1 

1 



we have a situation in which (it is easy to show that) there exists a 
solution for all B £ E 2 . Of course here for some choices of B the solution 
is not unique. 

2.5 Theorem 2 
Let P and Q denote real n X n matrices such that 

Va > £ I Pu I and q u > £ | g {i \ 

for all j = 1, 2, • • • , n. For j = 1, 2, • • • , n let /,-(•) denote a continuous 
monotone nondecreasing (but not necessarily differentiate) mapping of 
E 1 into itself, and let F(x) = fofa), / 2 (x 2 ), • • • , /„(x„)] tr for all x £ FT. 
Then for each R £ E n , there exists a unique x £ FT such that 

PF(x) + Qx = R, 

and, for any y £ E n , x is the limit of the sequence x co) , x (1) , • • • defined by 

y ln) = D P F(x ln) ) + D x (n> 

2/ (n+n + (P - D P )F(x M ) + (Q - D Q )x M = R 

for n ^ 0, in which D P and D Q are diagonal matrices whose diagonal 
elements coincide with those of P and Q, respectively. 
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2.6 Proof of Theorem 2 

Since the continuous mapping [D P F{-) + D Q ] of E n into E" possesses 
an inverse [D P F(-) + D ]~ l , the equation 

PF(x) + Qx = R 

possesses a unique solution x if and only if y = D P F(x) + D Q x is the 
unique solution of 

y + PF[(D P F(-) + D Q y l y\ + Q[(D P F(-) + D Q y l y] = R 

in which P = (P - D P ) and Q = (Q - D Q ). 

Therefore, by Banach's contraction-mapping fixed-point theorem, 
it suffices to show that with the metric p(y, z) = 2?_, I Vi — z * l> tne 
operator i/ denned by 

tffo) = PF[(D P F(-) + Do)" 1 !/] + 0[(DpF(.) + D,)-**] 

for all y E E", is a contraction mapping of E n into itself. We show this 
as follows. Let y G E n and z £ E"\ Using the fact that 

« - <W(<W*(0 + ^,) _1 «] + <W*[(<W/(0 + do/) _1 «] 

for all real a and all j = 1,2, • • • , n, in which d Pi and d , is the ;' th 
diagonal element of D P and D Q , respectively, it is a simple matter to 
verify that for all j: 



>', 



d Qi + dpiTi 
and 



(!/,- - Zj), 



(<W*(0 + rfoy) _1 2/y - (<W*(0 + dg#) l *i - dq . + df , r> . (J// - *i) 

in which r,- = 1 if y { = «,- , and, if ?/,- ^ z, , 

/,[(f/ P ,f,(-) + rf q ,-)"'.v,] - /,-[(<*/.,■/,■(•) + d Q f)"'gj] 
"' " (<W;(0 + do*)' 1 *! - (<W,(") + doi)"^ 

Thus 
fffo) - H(z) 



NONLINEAR EQUATIONS 109 

in which r t ^ 0. Therefore 

in which a Qj = £, v * I Q_n I and <r P/ = E.vi I />.» I- Since *qi < d Q> 
and o>, < r//,, for all j, there exists a positive constant < 1 such that 

/ op, + o>,-r/ \ 

max I- r-j 1 ^ p 

, W 0| + d P -Til 

for all r, ^ 0. □ 

2.7 Theorem 8 

If A G P and det A^O, # for each j = 1, 2, • • • , n: /,(•) is a con- 
tinuous mapping of E x into itself such that 

/,.(*,) = for all .r,- 

or 

/,(*,) > for all x, > c 

and 

ji(x,) < for all Xi < -c 

for some c ^ 0, then, with F(x) = [fi(x t ), f 2 (x 2 ), ■•■ , /„(a:„)] tr for all 
x E I'"', 

|| F(x) + Ax || -» co as ||*||-»». 

2.8 Proo/ o/ Theorem 3 
We note that 

|| F(a) + A a; || -* *> as || x || -» « 

if and only if 

|| A' l F{x) + x || — > « as || x || — * oo. 

With M - A~\ let 

il/F(x) +x = q. (14) 

Since A 6 Po , we have M 6 P •' Since il/ G Po , we have 1 for any 
y £ A" 1 and y 9* 6 

Vk(My) k ^ 
for some index A; such that y k 9^ 0. 
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Suppose that F(x) ?* 6. Then there exists an index k x such that 
UAx k J[MF(x)\ ht ^ 
with f kl (Xk,) 9* 0. Thus, using (14), 

UXx kl )[MF(x)] kl + f kl (x kt )x h , = f hl (x kl )q k , 
and 

Either x kl £ [— c, c] or not. If not, then fk t (x kt )x kl > and | x kl | ^ 
| q k , |. Therefore for some index k t , \ x kl | ^ 5 t = max (c, | g ftl |), whether 
or not F(x) = 6. 

Let M (kl> denote the matrix obtained from M by deleting the k Y 
row and column, and let M ikl) denote the k t column of M with the /b, 
entry removed. Similarly, let x {tl) , q lkl) and F ikl) (x ikl) ) denote the 
(n — l)-vectors obtained from x, q, and F(x), respectively, by removing 
the k t entry. Then 

M ikt) F {kl) (x lkl) ) +* (fa , = q {kl) - M lk J kl (x kl ). 

Since Af c * ,} £ P , we can repeat the argument given above. Thus 
there exists an index k 2 , different from k x , such that 

I x km | S S 2 = max (c, | q (kl , kl) |) 



in which 



I ?<*,.*.> I = max | [q (kl) - M l ''/t.OOL. 
i**,is«. 



and Z 2 is the index of the component of x (kl) that corresponds to the k% 
component of x. By continuing in this manner we can determine positive 
constants 5 X , 5 2 , • • • , 8 n depending only on q, F, M, and c such that, 
with 5 = max, {5,}, 

\xj \ ^ 8 for all ;' = 1, 2, • • • , n 

and each 5,- depends on g such that for any positive constant a, there 
exists a constant /3,-(a) with the property that 5, rjS jS«(a) provided that 
1 1 q 1 1 ^ a. Therefore for any a > there is a /3(a) such that 1 1 x \ | ^ /3(a) 
whenever \\q\\ ^ a, which implies that 1 1 ? 1 1 — * °° as 1 1 re 1 1 — > » . □ 

2.9 Theorem 4 

Let P and Q denote real n X n matrices with P strongly column-sum 
dominant. Suppose that there exists a real diagonal matrix D > such 
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that DP is strongly column-sum dominant and DQ is weakly column-sum 
dominant. Then jor all real diagonal matrices D, > and D 2 > 0, all 
eigenvalues of (PD X + QD 2 ) lie in the strict (that is, open) right-half 
plane. 

2.10 Proof of Theorem 4 

Since the strict right-half plane contains all of the eigenvalues of P, 
there exists choices of D t > and D 2 > such that every eigenvalue 
of {PD X + QD 2 ) lies in the strict right-half plane. Thus it suffices to 
show that (PA + QD 2 ) does not possess an eigenvalue on the boundary 
of the complex plane for all D, > and all D 2 > 0. In other words, it 
suffices to show that (with i = (—1)*) 

PD, + QD 2 + ital (15) 

is nonsingular for all D, > 0, all D 2 > 0, and all real constants u. 
Suppose that (15) is singular for some co and some Z>i > and some 
D 2 > 0. Then (DPD, + DQD 2 + iuD) is singular. But DPD, is 
strongly column-sum dominant and DQD 2 is weakly column-sum 
dominant. Thus M = (DPD t + DQD 2 ) is strongly column-sum domi- 
nant, and, since 

| ma + iudj \ > £_, | m if | 

for all j, in which d, is the j tb diagonal element of D, it follows that 
det (M + itaD) ^ 0, which is a contradiction. □ 

2.11 Definition 2 

With q and p nonnegative integers such that (p + q) > 0, let 3 de- 
note the set of all matrices M such that M = I, © M, M 9 © • • • 
© M, with 



AT* = 



and 



for all k = 1, 2, 





1 


* <""| 






1 


< <r J" < 


1 





< ■}« < 


1 



* As suggested, if q = 0, then .'/ = Mi © Af 2 
then jl/ = J . 



flf w while if p = 0, 
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2.12 Definition 3 

With q and p nonnegative integers such that (p + q) > 0, let 3(a) 
denote the set of all matrices M such that M = I a © M, © M 2 • • • 
©il/ p with 

1 -d? y 

-«J W 1 J 



,1/, = 



and 



for all k = 1,2, •■■ , p.* 
2.13 Theorem 5 



< 8 ( r k) ^ a ( r h) 
< fij w ^ a<*> 



Let f E 3, Zei H be a real matrix of order (2p + q), and suppose 
that M~ l H G P for all M G 3(a). Then 

(T + D 1 )~ l (H + D 2 ) G Po 
for all diagonal matrices D x ^ and Z) 2 ^ 0. 

2.14 Proo/ o/ Theorem 5 

Suppose that for some Dj ^ and D 2 ^ 

(T + DO" 1 ^ + £> 2 ) « ^o • 
Then there exists 2 a diagonal matrix D > such that 

(!T + DO-'CH + D 2 ) + Z) 
is singular. It follows that 

J7 + A + TD 
is singular, in which A = D 2 + Z^D. Since 

A + TD = M(A + D) 

in which il/ G 3(«)> it follows that 

# + M(A + D) 
is singular, and therefore that 



* Ah suggested, if q = 0, then M = il/, il/ 2 
then il/ = /,. 



M p , while if p = 0, 
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M^H + (A + D) 

is singular, which is a contradiction since 1 M' 1 H E P and ( A + D) is a 
diagonal matrix with positive diagonal elements. □ 

2.15 Theorem 6 

Let M~ l G E -Po for all M E 3, and let del G ^ 0. Let R be as defined 
in Section 1 4- Then for any TG3 

(T + D 1 y 1 [(I + GRY'G + D 2 ] E Po 
for all diagonal matrices D x ^ and D 2 ^ 0. 

2.16 Proof of Theorem 6 

Since det G ¥■ and M~ l G E Po for all M E 3, it follows (see the 
proof of Theorem 7 of Ref . 2) that 

M~\I + GRY l G E Po 

for all MG3. 

Suppose that for some TG3 and some D, ^ and _D 2 ^ 

(T + DO" 1 !^ + GRY'G + Z) 2 ] $ Po • 

Then, following the proof of Theorem 5, we would have 

det \M-\I + GRY'G + (A + D)\ = 

for some M E 3 and some diagonal matrix (A + D) with positive 
diagonal elements, which is a contradiction. □ 

III. ACKNOWLEDGMENT 

The writer is indebted to A. N. Willson, Jr. for carefully reading the 
draft. 

REFERENCES 

1. Sandberg, I. W., and Willson. A. N., Jr.. "Some Theorems on Properties of 

DC Equations of Nonlinear Networks." B.S.TJ., 48, No. 1 (January 1969), 
pp. 1-34. 

2. Sandberg, I. W., and Willson. A. N, Jr., "Some Network-Theoretic Properties 

of Non-Linear DC Transistor Networks," B.S.T.J.. 48, No. 5 (May-June 
1969), pp. 1293-1312. 

3. Varga, R. S., Matrix Iterative Analysis, Englewood Cliffs, New Jersey: Pren- 

tice-Hall, 1962, p. 23. 

4. Stern. T. E., Theory oj Nonlinear Networks and Systems, Reading, Mass.: 

Addison-Wcsley, 1965, pp. 42-43. 



114 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1970 

5. Palais, R. S., "Natural Operations on Differential Forms," Trans. Amer. Math. 

Soc, 92, No. 1 (1959), pp. 125-141. 

6. Holzmann, C. A., and Liu, R.. "On the Dynamical Equations of Nonlinear 

Networks with n-Coupled Elements," Proc. Third Ann. Allerton Conf. on 
Circuit and System Theory, U. of Illinois, 1965, pp. 536-545. 

7. Goldstein, A. A., Constructive Real Analysis, New York: Harper & Row, 

1967, pp. 41-15. 

8. Sandberg, I. W., "Some Theorems on the Dynamic Response of Nonlinear 

Transistor Networks," B.S.T.J., 48, No. 1 (January 1969), pp. 35-54. 

9. Hamming, R. W., Numerical Methods for Scientists and Engineers, New 

York: McGraw-Hill, 1962. 

10. Ralston, A. A., A First Course in Numerical Analysis, New York: McGraw- 

Hill, 1965. 

11. Hachtel, G. D., and Rohrer, R. A., "Techniques for the Optimal Design and 

Synthesis of Switching Circuits," Proc. of the IEEE, 55, No. 11 (Novem- 
ber 1967), pp. 1864-1876. 

12. Sandberg, I. W., and Shichman, H., "Numerical Integration of Systems of 

Stiff Nonlinear Differential Equations," B.S.T.J., 47, No. 4 (April 1968), pp. 
511-527. 

13. Calahan, D. A., "Efficient Numerical Analysis of Non-Linear Circuits," Proc. 

Sixth Ann. Allerton Conf. on Circuit and System Theory, U. of Illinois, 

1968, pp. 321-331. 



